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ABSTRACT 

Galactic collisions are normally modeled in a CDM model by assuming the 
DM consists of a small number of very massive objects. This note shows that 
the behaviour of a CDM halo during collisions depends critically on the mass 
of the particles that make it up, and in particular, all halo particles below a 
certain characteristic mass are likely to be lost. 

Subject headings: galaxies: dark matter, merging galaxies 
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1. Introduction 

There is now a sizeable literature on collisions between galaxies, and computer models 
show a remarkable resemblance to observations of interacting galaxies (see, e.g. [Barnes 



1982|) . Typically these involve the numerical solution of N-body collisionless Boltzmann 
equations, where N ~ 10 4 . This implies that the effective mass of a dark-matter (DM) body 
is 1O 6 M and one may wonder whether the results are independent of this assumption. 
In particular, it would be interesting to know whether the behaviour of cold dark matter 
(CDM) halo particles in a collision is effectively the same whether they are massive 
neutrinos or MACHOs. 

In this note we point out that CDM haloes may behave in radically different ways 
depending on the mass of the particles. Specifically haloes composed of particles lighter 
than a certain mass mo (mo ~ 10 16 kg) will be stripped off a galaxy during a collision, while 
heavier ones will remain with the galaxy and would presumably become part of a merged 
galaxy. We derive a specialized form of the Chandrasekhar dynamical friction formula 
which is the gravitational analogue of the Bethe-Bloch formula, and then apply this to a 
highly idealized galaxy collision. Finally we offer some speculations on how this will affect 
galaxy populations. 



2. Dynamical Friction 

Consider a particle of mass fi, which is incident on a gravitationally interacting gas 
consisting of particles of mass m and number density n. The dynamical friction formula 
( Chandrasehkar 1943] , can be written 



dt 
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We take the gas to be at rest, and the incident particle to be moving with a velocity v 
and simplify to the case m = p. None of these assumptions are critical to the argument. It 
is then easy to transform (1) to 

dE 8nG 2 p 2 p 



dx v 



2 



In (A) (2) 



where p is the mass density of the gas. The argument of the log is very large, so the log 
can be replaced by a dimensionless constant Cq. The resulting equation can be derived from 
the Bethe-Bloch (see, e.g. [Rossi 19521) equation by a judicious change of variables. Note 



the occurrence of the p 2 term, which is crucial to the following argument: qualitatively one 
can see that the cross-section is proportional to p and the mean momentum transfer is also 
proportional to p. It is also somewhat counter-intuitive that the force increases as \ at 
low velocities. To see how this affects a galactic halo, consider the collision process shown 
below 

EDITOR: PLACE FIGURE ?? HERE. 

. The galaxies are assumed to have a spherically symmetric halo of particles of mass p and 
density p(r), and to be moving with a relative speed of Voq. It is then trivial to integrate 
(2) to give 

v A = v A 00 - 167rC G 2 p J p{r)ds (3) 

where the integral is along the path followed by the particle. To proceed further we 
need to have a definite model for the density distribution in the galaxy. The model must 
give flat rotation curves out beyond the edge of the visible galaxy, and yet have a finite 
mass. For simplicity we choose 



^1 - e ( r °) 



M(r) = M 00 [l-e (4) 
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which gives a density function 

P(r)=p — ,y = — (5) 
y 2 r 

(This is, of course, not a good approximation for the core of a galaxy, but fortunately 
this is not relevant to the following argument). If the rotational speed of the galaxy is vo 
we get 

' P0 = 4^I (6) 

This allows (4) to be integrated giving 

v 4 = vi - 32nG 2 ^r C I (b ) , / (b ) = f t^— dy » ^ (7) 

i y\]y 2 - b 2 b ° 

where b = - is the scaled impact parameter and the approximation for the integral is 
good to about 50% for the range .01 < b < 2. By inspection of one can see that a particle 
will only be stopped in a galactic collision if its mass was larger than the characteristic 
mass [/,o given by 

Ho = V ~ r (8) 

2(27rG 2 p C /(6o)) 3 



3. An idealized Galaxy Collision 

. To give some idea of the value of /i , we can take plausible values of r = 50 kPc 
and the orbital velocity vo = 200kms~ 1 and further assumes that the galaxies have a 
velocity = lOO&ms -1 when they collide. This gives rise to = I0 42 kg ~ 10 12 M Q , 
p 10~ 22 kgm- 3 and C ~ 300 (for /i = 10~ 27 kg) to 30 (/i = 10 30 kg). This gives rise to a 
characteristic mass /io ~ 10 16 kg. Loosely speaking, any particle with a mass of less than 
/io would not be stopped in the collision. Note that this argument does not apply to gas 
or dust, which are affected by conventional friction. Hence we may expect that baryonic 
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matter would behave in the way predicted by the N-body simulations, and in particular 
MACHOs would not be lost to the merged galaxy. Possibly asteroids which are sufficiently 
massive that the conventional frictional force on them would be relatively small and which 
were not gravitationally bound to a solar system would also be lost. 

However, particles which interact only via gravitational interactions only, which means 
neutrinos or WIMPs (we use WIMPs generically in what follows) would behave very 
differently. WIMP haloes would be stripped off during a collision of this kind, with a time 
scale of — ~ 10 8 yrs. To see what would happen to the stars in such a collision, consider 



a collision between two identical spiral galaxies, which occurs in such a way that the cores 
coalesce immediately. Assuming that the luminous matter obeys a similar mass distribution 
function 



so that the galaxy is entirely luminous matter at the centre, and T\ ~ 10 kPc is the 
radius of the luminous matter. In this case, stars at the centre of the galaxies would be 
essentially unaffected by the loss of the halo, stars out to a critical radius R ~ 1.59ri would 
be forced into increasingly elliptical orbits and any stars beyond R would no longer be 
bound. 

The most serious weakness in this model is the assumption that the target gas is at 
rest. Since obviously nothing is known about the velocity of WIMP haloes, it is difficult 
to correct this. One might argue that the WIMPs form a gas with a BB temperature of 
(say) 2 K, which would correspond to velocities in the range of a few ms^ 1 , in which case 
the assumption is obviously valid. Alternatively the WIMP halo could co-rotate with the 
central galaxy, or (more plausibly) form a gravitationally self-bound system, in which case 
the virial theorem gives 




(9) 
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GM 00 ln(2) 



) 



(10) 



The latter two assumptions both imply speeds of 200 kms~ l whichare comparable to the 
collision speeds. This would require much more detailed modeling than is possible here, 
but one would expect that some fraction of the halo with a small velocity in the final rest 
system of the merged galaxy would be retained, while the majority would be stripped. 

The idea that haloes are stripped in collisions is in qualitative agreement with two 
results. Firstly |Dubinski et al 1996| show that tidal tails form during galaxy mergers appear 
to be considerably longer than would be predicted by N-body simulations. They note 
that the best estimates of dark matter derived from tidal tails in merging galaxies such as 
NGC 4038/39 and NGC 7252 imply a value about 1/3 of that derived from spiral galaxies. 
This could be understood if much of the halo material is lost after the collision. If, as is 
generally believed, elliptic galaxies are formed by collisions of spirals, they should be almost 
devoid of haloes. Galaxies in a cluster would have lost much of their DM haloes during 
repeated collisions, and the DM would be bound to the cluster as a whole. This is, of 
course, consistent with observations of (e.g.) the Coma cluster. 

One could regard this as strong circumstantial evidence that the main component of 
haloes are WIMPs. Obviously if this can be confirmed by more detailed models, it would 
be of considerable importance. 



4. 



Conclusion 



This calculation suggests that the behaviour of the DM haloes during a galactic 
collision depends critically on the mass of the particles that make it up. If the DM objects 
are relatively light (both neutrinos and WIMPs are light in this context), then the halo 
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will be lost totally, while heavy haloes would be retained as predicted by the N-body 
simulations. In other words, this process provides a way to differentiate between a particle 
type solution and a MACHO type solution to the DM problem. 

This work was suggested by a talk given by Jayanne English, and I am grateful to 
her for a useful correspondence. Conversations with Bob Carnegie and Steve Godfrey are 
acknowledged. This work was partly supported by NSERC. 
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